function meanerr(time,~,R,~,y0,Y)

load('orthbase');

 exactmean=zeros(1,100);
% time(1)=10^-15;

% for ii=1:100
%     exactmean(ii)=quadl(@(s) R(s),0,time(ii));
% end
%exactmean=y0(1)+exactmean-1/4*time.^2;
exactmean=-2*time;
hold on
plot(time,exactmean);
plot(time,Y,'r-');
title('Mean')
legend('exact mean','wce mean')
